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LONG-TERM  GOALS 

The  Navy  has  a  requirement  to  rapidly  and  covertly  characterize  the  coastal  environment  in  support  of 
Joint  Strike  Initiatives.  Over  the  past  16  years  we  have  demonstrated  that  spacebome  hyperspectral 
remote  sensing  is  the  best  approach  to  covertly  acquire  data  on  shallow  water  bathymetry,  bottom 
types,  hazards  to  navigation,  water  clarity  and  beach  and  shore  trafficability  to  meet  those 
requirements.  The  long  term  goal  of  this  work  is  to  put  a  hyperspectral  imager  capable  of  making  the 
appropriate  measurements  in  space  to  demonstrate  this  capability. 

OBJECTIVES 

The  objective  of  this  work  is  to  put  a  hyperspectral  imager  in  space  to  demonstrate  the  ability  to 
covertly  acquire  data  on  shallow  water  bathymetry,  bottom  types,  hazards  to  navigation,  water  clarity 
and  beach  and  shore  trafficability.  The  proposed  activities  are  designed  to  take  advantage  of  the 
Hyperspectral  Imager  for  the  Coastal  Ocean  (HICO)  currently  flying  on  the  International  Space  Station 
(ISS).  As  HICO  Project  Scientist  I  work  to  enhance  community  awareness  of  the  need  for  and  utility  of 
hyperspectral  imaging  of  the  coastal  ocean.  Our  work  includes  advancing  methods  of  on-orbit 
calibration  and  product  validation,  and  processing  and  analyzing  hyperspectral  data  of  the  coastal 
ocean.  In  particular  in  this  study  we  are  using  HICO  data  to  characterize  the  properties  of  river  systems 
and  coastal  waters.  Conventional  ocean  color  sensors  have  1  km  pixels  and  a  few  spectral  channels; 
these  have  not  proven  adequate  to  resolve  the  complexity  of  river  systems.  HICO  was  designed  to 
sample  the  coastal  ocean  and  has  95  m  GSD  and  88  spectral  channels  suitable  to  resolve  river  systems 
and  other  coastal  features.  To  make  the  best  use  of  this  data  we  are  developing  algorithms  and 
approaches  for  resolving  the  constituents  of  river  plumes,  harmful  algal  blooms  and  other  complex 
systems. 

APPROACH 

Two  tasks  are  planned  for  the  two  year  period  of  this  grant  (FY2010-201 1): 

1 .  Develop,  test  and  evaluate  algorithms  for  deriving  optical  properties,  chlorophyll,  suspended 

sediments  and  CDOM  for  coastal  systems  including  river  systems.  We  will  work  with  MERIS  and 
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HICO  data  and  use  the  Columbia  River  system  and  adjacent  coastal  waters  as  a  test  area  for  this 
work. 

2.  As  soon  as  it  is  available  we  will  collect  HICO  data  for  the  Yangtze  River  and  adjacent  coastal 
regions  in  China.  Once  we  have  developed  the  algorithms  and  approaches  that  work  for  the 
Columbia  River  we  will  test  and  evaluate  those  algorithms  for  the  Yangtze  and  other  river  systems. 
The  goal  is  to  validate  our  algorithms  and  to  further  our  understanding  of  this  important  river  and 
the  East  China  Sea  which  are  rapidly  changing  due  to  the  development  of  the  Three  Gorges  Dam 
and  continued  urbanization  of  China. 

WORK  COMPLETED 

Our  work  is  focused  on  the  coastal  ocean  and  a  major  issue  for  the  coastal  ocean  is  that  the  standard 
case  1  algorithms  used  to  calculate  chlorophyll  and  other  water  properties  assume  that  phytoplankton 
with  an  associated  level  of  Colored  Dissolved  Organic  Matter  (CDOM)  and  water  itself  are  the  only 
optically  active  components.  In  coastal  waters  high  levels  of  CDOM  from  rivers  and  coastal  runoff, 
large  phytoplankton  blooms,  sediments  from  rivers,  or  resuspension  from  the  bottom  are  all  significant 
optical  components  that  need  to  be  considered  as  part  of  the  optical  signature  (Davis,  et  al.,  2007).  For 
example  the  standard  MODIS  and  MERIS  products  give  false  high  chlorophyll  values  for  the 
Columbia  River  Plume  (Fig.  1).  The  MERIS  neural  network  (algal  2)  algorithms  are  designed  for 
European  coastal  waters  and  do  a  better  job  of  separating  chlorophyll  and  suspended  sediments.  We 
are  working  to  modify  them  for  Oregon  coastal  waters  and  eventually  for  use  with  HICO  data.  Note 
also  that  the  1  km  MODIS  data  does  not  adequately  sample  the  Columbia  River  mouth  including  the 
mixing  zone  that  is  order  50  km  inland  from  the  coast.  MERIS  300  m  data  does  a  better  job  of  imaging 
the  estuary  but  HICO  has  95  m  GSD  and  full  spectral  data  (88  channels  covering  the  400  -  900  nm 
spectral  region)  for  this  example  river  system. 


Figure  1.  Satellite  ocean  color  images  of  the  Columbia  River  estuary  and  plume  on  September  10, 
2009.  Left,  is  the  MODIS  standard  1000  m  chlorophyll  product.  This  algorithm  gives  false  high 
values  for  the  river  system  (red  in  this  color  scale)  due  to  the  high  suspended  sediments  in  the  river 
water.  The  large  pixels  do  not  image  the  river  mouth  effectively.  Center  is  the  MERIS  300  m 
standard  chlorophyll  product  (algal  1).  The  river  mouth  and  near  shore  plume  are  black  indicating 
the  algorithm  does  not  give  a  valid  product  for  these  waters  and  the  results  are  masked  out.  Right  is 
the  MERIS  neural  network  coastal  chlorophyll  product  (algal  2)  which  shows  reasonable 
chlorophyll  values  for  the  river  mouth  and  nearshore  plume.  The  300  m  MERIS  pixels  do  a  much 

better  job  of  imaging  the  river  mouth. 
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The  Hyperspectral  Imager  for  the  Coastal  Ocean  (HICO;  Lucke  et  al.  2011;  Corson  and  Davis,  201 1)  is 
an  imaging  spectrometer  based  on  the  PHILLS  airborne  imaging  spectrometers  (Davis  et  al.  2002). 
HICO  is  the  first  spacebome  imaging  spectrometer  designed  to  sample  the  coastal  ocean.  HICO 
samples  selected  coastal  regions  at  95  m  with  full  spectral  coverage  (400  to  900  nm  sampled  at  5.7  nm) 
and  a  high  signal-to-noise  ratio  to  resolve  the  complexity  of  the  coastal  ocean.  HICO  is  sponsored  by 
the  Office  of  Naval  Research  as  an  Innovative  Naval  Prototype  (INP),  to  demonstrate  coastal  products 
including  water  clarity,  bottom  types,  bathymetry  and  on-shore  vegetation  maps.  As  an  INP,  HICO  also 
demonstrates  innovative  ways  to  reduce  the  cost  and  schedule  of  this  space  mission  by  80%  by 
adapting  proven  PHILLS  aircraft  imager  architecture  and  using  Commercial  Off-The-Shelf  (COTS) 
components  where  possible. 

The  HICO  program  was  initiated  in  February  2006.  In  January  2007  HICO  was  selected  to  fly  on  the 
Japanese  Experiment  Module  Exposed  Facility  (JEM-EF)  on  the  International  Space  Station. 
Construction  began  following  the  Critical  Design  Review  on  November  15,  2007.  HICO  was 
completed  in  July  2008  and  it  was  integrated  into  the  HICO  and  RAIDS  Experimental  Payload  (HREP) 
in  August  2008.  HICO  is  integrated  into  HREP  and  flown  with  support  and  direction  from  DOD’s 
Space  Test  Program.  HREP  was  launched  on  the  Japanese  H-2  Transfer  Vehicle  (HTV)  September  10, 
2009.  The  HTV  rendezvoused  with  the  ISS  on  September  17,  2009.  HICO  was  installed  on 
September  23  and  collected  its  first  images  on  September  24,  2009.  To  date  over  4000  HICO  scenes 
have  been  collected,  procesed  to  level  lb  and  archived.  The  researech  community  can  access  HICO 
data  throught  the  OSU  HICO  website:  http://hico.coas.oregonstate.edu. 

For  the  past  two  years  we  have  been  collecting  HICO™  data  for  the  Columbia  (Fig.  2)  and  Yangtze 
Rivers  (Fig.  3).  There  are  many  constraints  on  data  collection  with  this  demonstration  instrument 
including  being  limited  to  one  image  per  orbit  for  data  transmission,  gaps  due  to  the  ISS  orbit  and 
operations  and  clouds.  However,  we  now  have  a  number  of  good  images  for  both  rivers  and  are 
continuing  to  collect  imagery  at  both  locations. 

Yangtze  River  in  China  is  a  major  source  of  sediments  and  nutrients  to  the  China  Sea  and  Straits  of 
Taiwan.  In  a  comparison  of  chlorophyll  products  from  MODIS  and  HICO™  (Fig.  3)  we  see  the  big 
picture  in  the  MODIS  data,  but  far  more  detail  in  the  HICO™  data.  Also,  note  that  the  MODIS 
algorithms  fail  over  waters  with  high  sediments,  but  the  HICO™  data  is  not  saturated  and  the  HICO™ 
algorithm  returns  useful  data  even  in  these  high  sediment  waters.  The  at  sensor  radiances  are  very 
close  and  were  the  MODIS  data  is  not  saturated  the  match  of  chlorophyll  values  is  very  good. 
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Figure  2.  Psudocolor  image  of  the  Columbia  River  mouth  and  adjacent  coastal  waters  made  with  3 
channels  of  HICO  data.  The  90  m  HICO  data  shows  many  more  details  including  shallow  bars  in 
the  river,  the  bridge  from  Astoria  Oregon  across  to  Washington,  breakwaters  and  the  complex 
mixing  patterns  as  the  plume  moves  offshore.  HICO™  is  a 
trade  mark  of  the  US  Naval  Research  Laboratory. 


Figure  3.  Nearly  coincident  MODIS  and  HICOTM  images  of  the  Yangtze  River,  China  taken  on 
January  18,  2010.  Left,  MODIS  image  ( 0500  GMT)  of  Chlorophyll-a  Concentration  (mg/m3) 
standard  product  from  GSFC.  White  areas  are  regions  where  sensor  channels  saturated  and  the 
algorithm  failed  to  produce  a  product.  The  box  indicates  the  location  of  the  HICO  image  relative  to 
the  MODIS  image.  Middle,  HICO  image  ( 0440  GMT)  of  Chlorophyll-a  Concentration  (mg/m 3) 
from  HICO  data  using  ATREM  atmospheric  correction  and  a  standard  chlorophyll  algorithm.  Left, 
comparison  of  the  at  sensor  radiances  from  HICO  (red)  and  MODIS  (blue)  for  location  X  in  HICO 
image  showing  the  accurate  on-orbit  calibration  for  HICO. 

( Unpublished  results  by  R-R  Li  and  B-C  Gao.) 


4 


To  validate  MERIS  and  HICO  data  products  for  the  Oregon  Coast  we  have  been  collecting  profiles  of 
downwelling  irradiance  and  upwelling  radiance  using  a  Satlantic  HyperPRO  (Fig.  4).  The  HyperPRO 
is  a  free  falling  optical  profiling  system  that  collects  profiles  of  spectral  Lu  and  Ed  and  chi 
fluorescence,  backscatter,  T  and  salinity.  The  system  is  calibrated  by  Satlantic  and  we  use  the  Satlantic 
software  for  processing  including  all  of  the  latest  corrections  based  on  NIST  calibrations. 


Figure  4.  Collecting  HyperPRO  data  and  an  example  data  set  from  the  MILOCO  cruise  off  the 
Oregon  coast  taken  June  4,  2009.  (Left  side  of  the  figure  shows  a  research  associate  lowering  the 
HyperPRO  instrument  over  the  side  of  a  small  research  vessel.  The  center  panel  shows  the  depth 
profile  of  the  spectra  of  downwelling  irradiance.  The  right  panel  shows  the  depth  profile  of  the 
spectra  of  the  upwelling  radiance  collected  with  the  HyperPRO  instrument.) 


This  system  produces  high  quality  measurements  of  spectral  remote  sensing  reflectance  (Rrs)  for  direct 
comparison  to  the  HICO  data  after  atmospheric  correction.  The  HyperPRO  data  together  with  other 
data  collected  on  each  station  including  HPLC  pigments,  productivity,  CDOM,  suspended  sediments 
are  placed  in  the  MILOCO  data  base  with  web  access. 

A  key  focus  of  our  effort  is  to  differentiate  the  river  plume,  Harmful  Algal  Bloom  (HAB),  or  other  key 
spectral  feature  from  the  background  signal  and  to  find  a  rapid  way  of  processing  the  data  to  produce  a 
river  plume  product.  While  standard  methods  exist  for  atmospheric  correction  of  land  (Gao,  et  al. 
1993)  and  open  ocean  (Gordon  and  Wang,  1994)  data  atmospheric  correction  of  coastal  ocean  data 
remains  problematic.  With  that  in  mind  our  method  uses  techniques  that  do  not  need  atmospheric 
correction,  but  start  directly  with  the  calibrated  at-sensor  radiances. 

The  basic  idea  for  the  indicator  methods  is  to  first  estimate  what  a  typical  'dark  pixel’  spectrum  for  a 
region  of  interest,  and  second  (after  subtracting  this  'typical  dark  pixel’)  to  examine  the  sensitivity  of 
the  residual  spectrum,  including  spectral  derivatives,  to  target  products  —  sediment,  chlorophyll,  or 
sampled  pixels  know  to  contain  pigments  of  interest,  such  as  phycocyanin  commonly  found  in 
cyanobacteria  associated  with  HABs.  That  is,  the  indicators  are  meant  to  be  'fingerprints’  for  the  target 
product  (Bustillos-Guzman  et  al.  2004). 
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The  starting  point  is  converting  the  at-sensor  radiance  data  to  a  normalized  reflectance  by  converting  to 
the  apparent  at-sensor  reflectance.  This  is  computed  using  the  Naval  Research  Lab’s  tafkaa_6s  code 
(Gao,  et  al.  2000)  according  to  the  formula: 


Pobs  /j.qE0  (1) 

where  Lt  is  the  observed  radiance,  Eq  is  the  solar  irradiance,  and  po  is  the  cosine  of  the  solar  zenith 
angle.  Next  we  create  a  reference  spectrum  for  each  image.  More  specifically,  we  try  to  estimate  a  dark 
ocean  pixel  from  looking  at  the  sensor  radiance, 

r  *  obs(l )  =  r  a(l )  + 1  r  w(l )  (2) 

where  r  a  accounts  for  atmospheric  and  sea-surface  reflection,  r  w  picks  up  contributions  below  the 
water  surface,  and  t  denotes  the  transmittance  from  the  water  surface  to  sensor.  In  developing 
‘indicator  maps,’  we  will  consider  a  slightly  different  decomposition.  Namely,  we  start  by  imagining,  at 
every  pixel,  the  at-sensor  radiance  we  would  see  if  the  water  was  clear;  a  so-called  dark  water  pixel. 
This  is  an  ‘idealized’  quantity,  but  the  difference  between  this  ideal  radiance  signal,  and  the  observed 
signal  at  sensor,  provides  information  about  what  is  in  the  water,  which  is  relatively  independent  of 
what  is  in  the  atmosphere.  Thus  it  is  a  good  starting  point  for  creating  an  ‘indicator’  function  for  what 
is  in  the  water. 

In  our  approach  a  dark  water  pixel  will  be  defined  ‘empirically,’  based  on  an  image,  or  collection  of 
images,  of  a  region  of  interest,  and  choosing  pixels,  or  patches  of  pixels,  which  determine  ‘dark  water,’ 
or  water  free  of  river  plume  materials,  or  other  features  of  interest,  for  that  region.  So  in  practice,  the 
dark  water  pixel  could  contain  some  background  material  in  the  water  column  that  we  will  consider  as 
‘clear  water,’  or  a  dark  water  pixel,  which  is  typical  for  the  region.  We  call  this  type  of  pixel,  which 
contains  some  below  water  signal,  a  regional  dark  pixel,  but  in  the  following  discussion  we  will  refer  to 
is  as  a  dark  pixel.  It  is  the  background  signal  from  which  we  start  any  further  signal  processing. 


Mathematically  we  can  write  this  decomposition  as: 

(3) 

(4) 


r*obs(l  )  =  ra(l  )  +  t  (r d(l  )  +  rb(l )) 
=  r  M(l )  +  r  i(I  ) 


where  rM(l  )  is  the  modeled  spectrum,  and  r  i(l  )  is  the  ‘indicator’  spectrum,  it  is  simply  the  residual 
between  the  at  sensor  radiance  minus  any  modeling  we  do  for  the  spectrum  (Tufillaro  and  Davis  2010). 
If  r  M  is  simply  chosen  as  a  dark  pixel  in  the  scene  than  it  is  what  is  normally  referred  to  as  ‘dark  pixel 
subtraction’,  a  very  simple  but  often  effective  scheme  for  atmospheric  correction.  The  terms  r  d(l  )  and 
r  b(l  )  are  called  the  ‘dark  water’  and  ‘bright  water’  contributions  to  r  w(l  )  respectively.  The  trick  to  this 
approach  is  finding  (a  typically  empirical)  model  of  dark  pixel  spectra  that  allow  us  to  create  a  data 
based  decomposition  of  r  *0bs(l  )  to  identify  one  or  more  dynamic  water  constituents  of  interest. 
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In  the  visible  spectrum  the  dominant  spectral  feature  is  Rayleigh  scattering  which  has  the  spectral  form 
of  a  negative  power  function.  To  create  an  empirical  model  to  fit  the  ‘dark  water’  pixels  we  take  a 
guess  at  the  following  functional  form: 


r  m(I  n)  -  (a  +  b(l  n- 1  i))’4 


(5) 


where  a  and  b  are  model  parameters  which  are  to  be  estimated  from  the  spectral  data  2„  and  r  M(l  n)  and 
2/  is  the  first  value  of  the  wavelength  in  data  set  being  modeled,  a  constant.  In  HICO  Lib  data  sets, 
each  at-sensor  pixel  has  88  radiance  values  between  400  nm  to  900  nm.  The  typical  scene  size  covers 
approximately  42  km  by  190  km,  and  500  by  2000  pixels.  To  ‘model’  the  data  we  limit  our  sets  to 
wavelengths  between  450  nm  to  900  nm,  so  in  our  data  sets  2/  =  450  nm  and  An  =  900  nm.  We  have 
applied  it  to  the  Columbia  River  and  Yangtze  River  HICO  images  using  it  to  differentiate  the  plume 
features.  Those  results  are  shown  in  the  results  section  below. 

The  next  step  is  to  take  the  radiances  (less  dark  pixels),  compute  its  derivatives,  and  then  create  a  map 
(typically  just  a  nonlinear  regression)  from  the  reflectance  data  (r  w,  d  r  w  /d2,  d  r  w  / d  2, ...)  to  our 
target  product,  either  pixels  in  the  scene  known  to  contain  the  product  of  interest  or  the  spectra  of  in- 
situ  above  water  radiances  for  the  site  (Tufillaro,  Davis,  and  Ortiz  201 1).  In  practice  we  often  compute 
the  'principal  components’  of  the  spectral  signal  first,  and  use  these  to  create  the  indicator  map 
following  the  procedure  recently  described  by  Oritz  et  al.  (Oritz  et.  al.  2011).  The  principal  components 
are  spectral  decomposition  of  the  signal  which  breaks  the  signal  into  its  empirical  modes  producing  the 
greatest  variance.  In  this  way  we  can  both  reduce  the  dimension  of  the  map,  as  well  as  pick  out  the 
parts  of  the  signal  which  are  best  correlated  to  changes  across  the  scene. 

The  use  of  derivative  data,  in  addition  to  the  original  spectral  radiances  is  helpful  since  it  can  highlight 
significant  features.  This  observation  has  been  used  often  in  the  past  for  identification  of  laboratory 
spectra,  and  it  is  possible  to  directly  adapt  laboratory  identification  methods,  such  as  derivative 
spectroscopy,  to  remote  sensing  product  estimation.  Derivate  spectroscopy  has  previously  been  used  to 
identify  optimal  spectral  channels  for  the  design  of  multi-spectral  instruments  (Lee  et  al.  2007),  and  to 
estimate  coastal  bottom  types  from  aerial  hyperspectral  data  (Louchard  et  al.  2002).  Here  we  illustrate 
the  use  of  methods  adapted  from  derivative  spectroscopy  for  product  generation  from  HICO  data. 
Derivate  spectroscopy  methods  are  useful  in  untangling  spectral  components  when  the  underlying 
scattering  or  absorptive  features  have  significantly  different  half  widths.  Consider  the  sketch  of  a 
simple  at-sensor  spectrum  shown  in  Fig.  5.  The  broadest  feature  is  the  Rayleigh  scattering  which  is  the 
monotonically  decreasing  across  the  spectrum  from  blue  to  red.  Taking  the  first  derivative  will  de- 
emphasize  this  broad  scattering  signature  by  essentially  subtracting  a  baseline  from  initial  spectrum. 
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Figure  5.  Sketch  of  a  typical  HICO  spectrum  showing  how  derivative  spectroscopy  can  amplify 
narrow-band  features  in  spectrum.  The  4th  derivative  of  the  spectrum  in  red  shows  how  the  fourth 
derivative  is  much  more  sensitive  to  narrow  spectral  features,  such  as  chlorophyll  fluorescence. 

Features  seen  include  a  signal  in  the  yellow  part  of  the  spectrum  with  a  relatively  board  bandwidth,  and 
a  signal  in  the  red  (685  nm)  for  chlorophyll  fluorescence  with  a  narrower  bandwidth.  The  detection  of 
these  underlying  absorptive  and  fluorescent  signatures  can  be  sharpened  by  the  use  of  derivative 
spectroscopy.  Specifically  peaks  with  narrower  half-widths  grow  more  quickly  with  the  order  of  the 
derivative.  This  is  illustrated  in  Fig.  6  where  we  show  that  the  derivative  helps  to  highlight  signatures 
sensitive  to  sediment  in  recent  Columbia  River  spectra.  Spectral  features  identified  in  this  way  can 
create  products  which  are  'regionally  tuned,’  and  built  on  historical  data  specific  to  a  coastal  area. 


(a)  Derivative  of  HICO  Spectra  (b)  HICO  Columbia  River  Swath  23  July  201 1 


wavelength 
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Figure  6.  Analysis  of  a  HICO  spectra  using  derivative  spectroscopy,  (a)  The  derivative  spectra 
calculated  at  the  points  indicated  in  the  HICO  image  (b)  taken  23  July  2011.  (c)  A  'regionally 
tuned ’  map  of  sediment  concentration  (bright  red  is  high  sediment,  dark  red  low) 

emanating  from  the  Columbia  River. 
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The  method  for  our  sediment  and  pigment  products  from  HICO  is  similar  to  that  recently  described  by 
Gitelson,  et  al.  (2011).  They  used  HICO  data  to  create  regionally  tuned  chlorophyll-a  product  for  the 
Azov  Sea,  where  they  pick  the  bands  for  the  product,  starting  with  all  the  HICO  wavelengths,  and  then 
optimizing  the  choice  of  bands,  for  a  band-ratio  algorithm,  based  on  in-situ  data.  Their  method  is  also 
informed  by  derivative  analysis  in  their  choice  of  the  final  bands  for  the  regression.  Our  method  differs 
in  that  we  build  our  optimization  using  all  the  HICO  wavelengths,  and  instead  of  testing  for  the 
sensitivity  to  specific  bands,  rather,  our  sensitivity  is  based  on  the  principal  components  which  are  a 
linear  combination  of  all  the  wavelengths.  Specifically,  the  method  we  describe  could  contain  details 
about  the  spectral  shape  -  such  as  spectral  bandwidths  -  which  might  not  be  utilized  in  the  method  of 
Gitelson  and  co-workers. 

Additionally  as  HICO  project  scientist,  I  am  funded  by  NRL  to  continue  to  work  with  the  engineers  and 
scientists  at  NRL  and  partner  institutions  to  prepare  for  the  processing  and  analysis  of  HICO  data.  The 
NRL  team  processes  the  HICO  data  to  standard  level  lb  (calibrated  at  sensor  radiances)  products.  At 
OSU  we  archive  and  distribute  HICO  data  for  academic  users  and  international  partners.  We  are  also 
processing  the  data  on  a  request  basis  providing  some  L3  products  with  geolocation  and  atmospheric 
correction,  as  well  as  collaborating  with  a  number  of  partners  to  use  HICO  data  as  part  of  their  on¬ 
going  programs. 

While  HICO  is  now  flying  on  the  ISS  I  continue  to  support  the  Naval  Research  Laboratory  in  the  effort 
to  fly  COIS  (Wilson  and  Davis,  1999)  on  a  spacecraft  of  opportunity  to  provide  higher  resolution 
hyperspectral  data  for  the  coastal  ocean.  Each  year  we  have  presented  COIS  to  the  Navy  and  DoD 
Space  Experiment  Review  Boards  (SERBs)  and  received  high  rankings.  While  HICO  data  is  95  m 
GSD  optimal  for  characterization  of  optically  deep  coastal  waters,  COIS  would  provide  the  30  m  GSD 
data  optimal  for  optically  shallow  waters  producing  products  like  bathymetry  and  bottom  types. 

RESULTS 

The  focus  of  this  effort  is  on  the  processing  and  analysis  of  HICO  data  to  demonstrate  the  utility  of 
hyperspectral  imaging  for  characterization  of  coastal  waters.  Our  efforts  are  particularly  focused  on 
river  systems  and  we  have  conducted  initial  analysis  of  HICO  data  for  the  Columbia  and  Yangtze 
Rivers.  During  this  first  year  we  have  processed  an  extensive  set  of  MERIS  data  for  the  Columbia 
River  and  Oregon  Coastal  Waters.  We  have  also  collected  an  extensive  set  of  in  situ  data  for 
validation  of  products  off  the  Northwest  Coast  and  are  currently  comparing  the  MERIS  products  with 
that  in  situ  data.  It  became  clear  early  that  both  full  resolution  (FR,  300  m)  MERIS  and  HICO  show 
great  potential  for  tracking  harmful  algal  blooms  (HAB’s)  Fig.  7  shows  the  use  of  MERIS  FR  (300m) 
data  to  track  an  intense  algal  bloom  during  September  2011. 

We  are  able  to  get  an  even  better  look  both  spatially  and  spectrally  using  HICO,  as  shown  in  Fig.  8. 

The  water  around  this  particular  bloom  is  filled  with  chlorophyll  according  to  standard  processing. 
However,  the  areas  indicated  by  red  west  of  the  Juan  de  Fuca  Strait  are  the  site  of  a  Pseudo-nitzchia 
bloom  (confirmed  by  in  situ  measurements).  Such  intense  surface  blooms  are  often  recognizable  as  a 
peak  that  is  measured  in  the  MERIS  709  nm  channel.  The  data  from  HICO  is  allowing  us  to  create  an 
even  finer  spectral  finger  print  of  the  diatoms;  in  this  case  Pseudo-nitzchia  which  often  produces 
HAB’s  in  the  Northwest  waters. 
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MERIS  RGB  201 1  09  04 


MERIS  Example  Spectrum  at  intense  surface  bloom 


Figure  7.  Tracking  an  intense  algal  bloom  down  off  the  northwest  coast  using  a  sequence  of five 
MERIS  FR  images  during  September  2011.  The  line  height  of  the  MERIS  Lib  709  nm  channel  is 
used  to  compute  the  Maximum  Chlorophyll  Index  (MCI)  to  map  the  algal  bloom. 
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Figure  8.  A  map  of  the  Pseudo-nitzchia  bloom  off  the  Northwest  coast  during  September  of  2011 
created  with  a  spectral  identification  fingerprint  from  HICO  Lib  data.  (HICO  image,  blow  up  of  the 
bloom  area,  and  comparison  of  spectra  of  the  surrounding  water  and  the  bloom 
water  with  a  spectral  peak  at  709  nm. 


Work  continues  on  processing  HICO  images  to  fine  tune  our  spectral  signature  methods  for  a  variety  of 
river  systems  including  the  Columbia,  Yangze,  Clutha  (NZ),  Han  (KR),  and  the  Amazon.  An  example 
of  the  where  we  have  applied  the  plume  indicator  approach  described  above  is  shown  in  Fig.  9  for 
HICO  images  of  the  Columbia  River.  More  recently,  we  began  work  with  Alex  Kurapov  (Physical 
Oceanography,  Oregon  State  University)  on  how  to  assimilated  ocean  color  data  of  the  Columbia  River 
plume  from  both  MERIS  and  HICO  into  Alex’s  groups  regional  ocean  model  for  forecasting  off  the 
Oregon  Coast  (Kurapov  et.  al  2011).  This  work  should  allow  us  to  make  short  term  forecasts  of  water 
clarity  and  turbidity  in  the  region  of  the  Columbia  River  mouth,  and  to  follow  the  fine  details 
(temporally  and  spatially)  of  the  wanderings  of  the  river  plume. 
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Figure  9.  Images  of  Columbia  River  with  indicator  map  highlighting  sediments:  (a)  LIB  HICO™ 
image  of  Columbia  River  mouth  19  March  2010.  (b)  An  indicator  for  sediments  in  the  Columbia 
River  based  on  the  dark  pixel  spectral  separation  for  the  image  presented  in  Fig.  3(a).  (c-f)  Typical 
spectra  and  sediment  indicator  maps  for  three  HICO  images  from  July  2010. 


Turning  to  the  next  example,  we  examine  the  Yangtze  River  using  HICO.  Fig.  10  shows  a  sequence  of 
clear  images  from  HICO  along  with  typical  plots  of  the  second  derivative  of  the  spectrum  with  respect 
to  wavelength.  The  wavelengths  around  620  are  sensitive  to  sediment  concentrations  (and  580  is  a 
chlorophyll  absorption  minimum),  and  MERIS  uses  a  620  nm  channel  to  estimate  sediment.  As  also 
pointed  out  by  Gitelson  et.  al.  (2011),  HICO  allows  us  to  tune  the  product  algorithm  wavelengths  to  the 
products  for  a  particular  region.  In  the  Yangtze  region,  which  can  exhibit  extremely  high  sediment 
concentrations,  the  sediment  maximum  occurs  closer  to  610  nm  rather  than  620nm,  and  our  principal 
component  method,  like  optimized  band  ratio  methods  (Gitelson  et.  al.,  1  201 1),  will  weigh  the 
sediment  product  algorithm  more  closely  to  the  sediment  maximum  when  trained  on  regional  data  (Fig. 
11). 
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Figure  10.  HICO  images  and  derivative  spectra  of  the  Yangtze  river  in  China.  Using  collections  of 
images  we  are  building  up  signatures  to  distinguish  the  constituents  of  the  water  column. 


Figure  11.  Upper,  HICO  RGB  image  and  sediment  product  map  for  Yangtze  River  in  China  on  6 
July  2010.  Lower,  the  second  derivative  of  the  spectrum  indicates  that  wavelenths  around  605  are 
much  more  sensitive  to  sediment  concentration  than  the  620  nm  MERIS  band  6.  Our  HICO 
product  algorthim  will  automatically  optimize  the  product  algorthim  to  weight  data  more  heavily  to 
around  605  nm  when  trained  on  Yangtze  regional  historical  data,  like  that  presented  in  Fig.  10. 
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HICO  has  been  on  the  ISS  for  two  year.  Over  4000  images  have  been  collected  for  locations  around  the 
world.  HICO  is  operating  as  planned  and  the  data  products  look  very  good.  We  continue  to  work  with 
the  NRL  team  on  the  on-orbit  calibration  and  validation  of  HICO  data.  As  Project  Scientist  it  is  my 
goal  to  work  with  the  scientific  community  to  make  the  best  possible  use  of  this  unique  data  set.  To 
that  end  we  operate  the  HICO  website  at  OSU  (http://hico.coas.oregonstate.edu )  which  serves  data  to 
interested  scientists  around  the  world.  We  currently  have  over  50  users,  a  number  of  them  have 
submitted  formal  proposals  and  their  projects  are  summarized  on  the  website.  The  website  also 
includes  publications  and  presentations  on  HICO,  an  archive  of  existing  data  and  a  tool  for  selecting 
sites  and  requesting  data  as  well  as  directions  for  working  with  the  data  and  products. 

IMPACT/APPLICATIONS 

The  long  term  goal  of  this  work  is  demonstrate  the  value  of  a  hyperspectral  imager  capable  of  making 
the  appropriate  measurements  in  space  to  demonstrate  the  capability  of  this  technology  for  the  rapid 
and  covert  characterization  of  the  coastal  ocean  to  support  naval  operations  around  the  world.  We  are 
using  data  from  HICO  on  the  ISS  to  demonstrate  that  capability.  The  work  completed  this  year  is 
another  incremental  step  towards  that  goal. 

RELATED  PROJECTS 

I  continue  to  collaborate  regularly  with  colleagues  at  the  NRL  Remote  Sensing  Division  (Code  7200; 
Mike  Corson  and  others)  and  the  NRL  Oceanography  Division  (Code  7300;  Bob  Amone  and  others) 
and  with  Zhong-Ping  Lee  now  at  U.  Mass.  Boston. 
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